Unified theory of quantum many-particle systems 

Yu-Liang Liu 

Center for Advanced Study, Tsinghua University, Beijing 100084, 
People 's Republic of China 

Abstract 

Using eigen-functional bosonization method, we study quantum many- 
particle systems, and show that the quantum many-particle problems end 
in to solve the differential equation of the phase fields which represent the 
particle correlation strength. Thus, the physical properties of these systems 
are completely determined by the differential equation of the phase fields. We 
mainly focus on the study of D-dimensional electron gas with/without trans- 
verse gauge fields, two-dimensional electron gas under an external magnetic 
field, D-dimensional boson systems, a D-dimensional Heisenberg model and a 
one-band Hubbard model on a square lattice, and give their exact (accurate 
for Heisenberg model) functional expressions of the ground state energy and 
action, and the eigen-functional wave functions of the fermions/bosons. With 
them, we can calculate a variety of correlation functions of the systems, such 
as single particle Green's functions and their ground state wave functions. In 
present theoretical framework, we can unifiably represent the Landau Fermi 
liquid, non-Fermi liquid (D > 2) and Tomonaga-Luttinger liquid. 
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I. Introduction 



Quantum many-particle systems are the main topics of the condensed matter physics, 
in which strongly correlated electron systems are the most interesting and hard problems, 
such as heavy fermion systems, high Tc cuprate superconductors, fractional quantum Hall 
effects, and some one-dimensional interacting fermion systems. In general, the quantum 
many-particle (fermion) systems can be divided into two categories, one is represented by 
the Landau Fermi liquid theory and may be called as weakly correlated systems, and 
another one is represented by non-Fermi liquid theory, such as the Tomonaga-Luttinger 
liquid theory [0-0] and marginal Fermi liquid theory |6|]7]], and may be called as strongly 
correlated fermion systems. Thus there explicitly exists a key parameter hidden in the 
quantum many-particle systems, which represents the fermion/boson correlation which is 
produced by the particle interaction. However, the traditional perturbation theories, such 
as Hartree-Fock approximation, mean field theory, random-phase approximation (RPA) and 
renormalization group theory, are starting from the bare Green's function to perturbatively 
treat the interaction (potential) terms of the particles by taking some low order Feynman 
diagrams, i.e., exactly treating the kinetic energy term and perturbatively treating the par- 
ticle interaction. In these methods, there does not exist a parameter used to represent the 
particle correlation which completely determines the physical properties of the strongly cor- 
related systems. Thus, the previous perturbation methods are successful in treating the 
weakly correlated systems, but using them to treat the strongly correlated systems, one can 
meet some serious problems, because in these systems it is difficult to find a suitable small 
quantity perturbation expansion parameter. 

It becomes clear that if one wants to unifiably represent the weakly and strongly cor- 
related systems, one must find a well-defined parameter which can represent the particle 
correlation strength. Due to the particle correlation is produced by the particle interac- 
tion, it is natural that we are starting from the interaction (potential) terms to exactly 
or perturbatively treat the kinetic energy term of the particles, i.e., we first exactly treat 
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the interaction terms, then with them we exactly or perturbatively treat the kinetic energy 
term. In this way, the basic ingredient of the systems is particle density field rather than 
particle field operators, like usual bosonization method in treating one-dimensional inter- 
acting fermion systems. In Ref. ||, we proposed an eigen-functional bosonization method, 
in which there naturally appears a key parameter called the phase field, and with it we can 
unifiably represent the weakly and strongly correlated electron gases, i.e., the Landau Fermi 
liquid theory and the non-Fermi liquid theory can be unified under our eigen-functional 
bosonization theory. 

According to the Hohenberg-Kohn theorem |J , the ground state energy of the quantum 
many-particle systems is uniquely determined by their ground state particle density. This 
theorem does not tell us how to construct the ground state energy by the particle density, 
but it clearly tells us that the ground state energy can be exactly represented by the particle 



density. However, in the Kohn-Sham scheme [10,11], we can only approximately obtain the 



expression of the ground state energy as a functional of the ground state particle density, 
because in this scheme the exchange correlation energy term is unclear. Only for small ex- 
change correlation energy can one use the Kohn-Sham scheme to obtain reliable results, i.e., 
one cannot use it to study the strongly correlated systems, where the exchange correlation 
energy is large . Thus it is very desirable to write out an exact expression of the ground 
state energy of the quantum many-particle systems as a functional of the particle density. 
In fact, in usual bosonization representation of one-dimensional interacting fermion systems 



I4IT3II , we learn how to construct the Hamiltonian by the fermion density operators, where 
we use the fermion density field to uniquely represent the kinetic energy and the fermion field 
operators. With the eigen-functional bosonization theory, we can treat any quantum many- 
particle system by using the fermion/boson density field to represent the kinetic energy and 



the eigen-functional wave functions of the fermions/bosons. In Ref. ||14fl , we demonstrated 
how to use the particle density field to exactly represent the ground state energy of the 
quantum many-particle systems. 

In this paper, we shall give a detail description of the eigen-functional bosonization 



method, and apply it to quantum many-particle systems in continuous coordinate space and 
some lattice models, such as the electron gases with or without transverse gauge fields, the 
two-dimensional electron gas under an external magnetic field, the interacting boson systems, 
the Heisenberg model and one-band Hubbard model. We clearly show that the problems of 
these systems end in to solve the differential equation of the phase fields which completely 
controls the physical properties of the systems, thus the quantum many-particle systems can 
be unifiably represented by the eigen-functional bosonization theory. In section II, we give 
the exact expressions of the ground state energy as a functional of the particle density field for 
different quantum many-particle systems. These expressions are universal, i.e., they are valid 
not only for weak fermion/boson interactions, but also for strong fermion/boson interactions. 
(1). we give the exact expression of the ground state energy as a functional of the electron 
density field for electron gases with and without transverse gauge fields, respectively; (2). 
we give the exact functional expression of the ground state energy for interacting boson 
systems. (3). we give the exact functional expression of the ground state energy of one-band 
Hubbard model on a square lattice, and the approximate functional expression of the ground 
state energy for the Heisenberg model. In section III, we focus on to calculate the action 
of the quantum many-particle systems. (1). we give the exact functional expression of the 
action for the electron gas with the transverse gauge fields, and the eigen-functional wave 
functions of the electrons. Under linearization approximation (only keeping the linear terms 
in solving the differential equation of the phase fields), we obtain its effective action which 
is the same as that obtained by usual random-phase approximation (RPA). (2). we give 
the exact functional expression of the action of the one-band Hubbard model on the square 
lattice, and the eigen-functional wave functions of the electrons. In section IV, we calculate 
the ground state wave function of the interacting boson systems, which is very similar to 
the correlated basis functions used in the study of the liquid 4 He. We also approximately 
calculate the ground state wave function of the two-dimensional electron gas under the 
external magnetic field, which has a little similarity with Laughlin's trial wave-functions 
T5| , |l6|l . In section V, by calculating the correlation function of the phase fields, we show 



that the imaginary part of the phase fields does represent the particle correlation strength. 
We also demonstrate that if the non-Fermi liquid behavior of the systems is derived from 
the electron Coulomb interactions, the correlation function of the phase fields has a log-type 
expression; If their non-Fermi liquid behavior is derived from the transverse gauge fields, the 
correlation function of the phase fields has a power-law form, t a , a < 1. Under our present 
theoretical framework, the Landau Fermi liquid, non-Fermi liquid and Tomonaga-Luttinger 
liquid can be unifiably represented. We give some discussions and our conclusion in section 



In this section, we shall give the exact expression of the ground state energy as 
a functional of the particle density field for a variety of quantum many-particle sys- 
tems. In general, a quantum many-particle system can be described by the Hamiltonian, 
i?T — Ho + H cx + i/gpin, where H presents the kinetic energy and the density-density in- 
teractions of the particles, H ex presents the interactions between the system and external 
fields, and H spin presents the spin-spin interactions among the particles, which is usually 
written out on lattice sites. In the following we consider the systems with these different 
particle interactions. 



For simplicity, we first consider an electron gas with the Hamiltonian (here omitting the 
spin label), 



where p = —ihV, p(x) = ip^(x)ip(x) is the electron density operator, and D the dimensions of 
the system. In usual perturbation theories, such as Hartree-Fock approximation, Random- 
Phase approximation (RPA) and renormalization group theory, one mainly focuses on how to 



VI. 



II. Exact functional expression of the ground state energy 



A. Electron gas 




(1) 
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perturbatively treat the density-density interaction (four-fermion interaction) term, and how 
to get more accurate results by considering some special Feynman diagrams. However, these 
methods do not always work for low- dimensional strongly correlated systems. According 
to the Hohenberg-Kohn theorem and the bosonization representation of one-dimensional 
fermion systems, we know that the ground state energy and/or the Hamiltonian can be 
uniquely represented by the ground state particle density, thus we can also study this system 
by perturbatively or exactly treating the kinetic energy term of the electrons. In order 
to calculate the ground state energy, we introduce a Lagrangian multiplier (or Hubbard- 
Stratonovich field) <j)(x) which takes p(x) = tp*(x)ij}(x) as a constraint condition, then the 
Hamiltonian ([]]) can be re-written as, 



where <p(x) and p(x) are independent boson fields, and the Hamiltonian H Q [(j),p\ only has 
the quadratic term of the electron field operator ip(x). It is clear that the kinetic energy 
and potential energy of the electrons are completely separated, and the electrons move 
in the "effective potential" <p(x). Only in this way (four-fermion interaction replaced by 
the "effective potential"), we can exactly treat the kinetic energy term. We would like to 
point out that after introducing the boson field 0(x), the Hilbert space of the system is 
enlarged, because the electron density field p(x) and the electron field operator now become 
independent of each other, and some unphysical electron density field p(x) 7^ i[)'(x)i/;(x) also 
appears in the Hamiltonian |2|. To only keep the physical electron density field p(x), we must 
add a constraint condition to the system (see below) SE[(f>, p]/5(j)(x) = 0. In analogy to the 
eigen-functional bosonization theory ||, here we solve the eigen- functional equation of the 
"propagator" operator (in correspondence with below M(x, t) which is a real propagator 
operator) M(x) = p 2 / (2m) — p + <p(x), 
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where *&k{x, [0]) are the eigen-functional wave functions of the electrons for a definite boson 
field <f>(x). Using the Helmann-Feynman theorem, we have the following expression of the 
eigen-values, 

E k [<j>]=e(k) + X k [<l>] 

£,[0] = J d^J d D x ( J ) (x)\y k (x,M)\ 2 (4) 

where e(k) = (hk) 2 /(2m) — \i. The eigen-functionals ^f k (x,[(f)]) have the following exact 
expression, 

m = ^* k '*e QkM (5) 

where A k is the normalization constant, and the phase fields Q k {x,^) satisfy usual (static) 
Eikonal equation, 

S + ^ ' p ) Qk(x ' Q + + mx) = El M (6) 

It is worthily noted that the eigen-functionals are composed of two parts, one represents 
the free electron, and another one represents the contributions from other electrons by the 
interaction potential v(x — y), therefore it is clear that the eigen-functionals are the eigen- 
functional wave functions of the interacting electrons corresponding to the definite boson 
field (x) . With them, we can construct the ground state wave-function, and calculate a 
variety of correlation functions of the system by taking a functional average over the boson 
field (x). 

The electron field operators ip^(x) and ip(x) can be represented as, respectively, 

^{x) = Y,n^W\)c\ 

k 

= 5>*(z, [0])c fc (7) 
k 

where c k (c k ) is the electron's creation ( annihilation) operator with momentum hk. With 
equations. ([5]) and (|7|), we can obtain the exact expression of ground state energy as a 
functional of the boson field 0(x) and the electron density p(x), 
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E[<f>,p] = f dPx < H [<j>,p] > 
J o 

= '£e(-E k [<f>])E k [<i>] - ! d D x<P(x)p(x) + i ! d D xd D yv(x - y)p(x)p(y) (8) 

k J J 

where the boson field is determined by the condition 6E[(/>, p]/S<j>(x) = 0, which leads 
to the equation, 



k 

where G k (x) = /„* d£e 2 G?(*.0/ / d D xe 2 ^ x ^\ and Q k (x,£) = +iQi(x,0- This 

equation shows that the boson field 4>(x) is the functional of the electron density field p(x). 
The chemical potential is determined by the constraint equation, 

J20(~E k m = N (10) 

k 

where N is the total electron number. With equation (f|), the ground state energy can be 
represented as another form, in which it is only the functional of the electron density field 
(4>(x) is the functional of p(x) determined by (H)), 

E g [p] =E + ~J d D xd D y (v(x - y)p(x)p(y) - 2<P{x)<P{y)^6{-E k [<P}) 5 -^^ (11) 

where Eq = 9(— E k [(f)\)e(k), and the self-energy can be written as a simple form, £&[</>] = 
J d D x(f>(x)G k (x). The equations @, (||), (|10D and (pT|) give the exact expression of the 
ground state energy as a functional of the electron density, we can numerically calculate 
it by taking SE g [p]/8p(x) = 0. These equations have more advantages than that in the 
Kohn-Sham scheme, because that, 1). the expression of the ground state energy as the 
functional of the electron density is exact, with it we can exactly determine the ground 
state energy and the ground state electron density; 2). with these equations, we can easily 
estimate the contributions from high order terms, and obtain enough accurate results we 
hoped for some special considerations; 3). Using these equations to calculate the ground 
state energy and the ground state electron density, we may need much less computer time 
than that in the Kohn-Sham scheme, because here we do not need to self-consistently solve 
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the eigen-equation of single electron. The ground state wave-function of the system can be 
obtained by the eigen-functionals, however, we cannot write it as a simple form, because the 
phase field Qk{x,C, = 1) is the function of the momentum k. 

The expression of the ground state energy (|11]) is valid not only for Landau Fermi liquid 
(weak correlation systems), but also for non- Fermi liquid (strongly correlated systems), 
because it is universal for weak and strong electron interactions. The phase field Qk(x,£) is 
a key parameter for unifiably representing the Landau Fermi liquid and non-Fermi liquid. 
Its imaginary part represents the electron correlation strength, and its real part (D > 2) 
only contributes to the ground state energy and the action of the systems, this can be clearly 
seen in the expression of the ground state energy (|TT|). 



B. Electron motion under an external magnetic field 

We now consider the electron motion under an external magnetic field B = (0,0, B), 
where the response of the system to the magnetic field cannot be written as a simple density 
form. It is well-known that for enough strong magnetic field, in low temperature limit a 
two-dimensional electron gas shows the fractional quantum Hall effects due to the Coulomb 
interaction of the electrons. Here we only give the exact expression of the ground state 
energy of this system, and do not compare it with that obtained by Laughlin's trial wave- 



functions [15], because if so it needs more detail numerical calculations. Under the magnetic 



field, the Hamiltonian (§) becomes, 

H[<f>, p] = ^{x) Q-(p + - c &) 2 -fi + (f>{x)j ij)(x) 

- (j){x)p{x) + ^ J d 2 yv{x - y)p{y)p{x) (12) 

where the gauge field A = (—yB/2,xB/2,0), and for simplicity we take D = 2. The 
eigen-functional equation corresponding to equation (§) reads, 

l-(p + e -Af-p + <P(x)\ * nl (x, [4>]) = E nl [^ nl {x, [<t>]) (13) 
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and the eigen-values are 

E n i\4>\ = e n + E n A 



J d£ j d 2 x<P(x)\^ nl (x,M)\ 2 (14) 



where e n = hu (n + 1/2) — /i, n = 0, 1, 2, loq = eB / (mc) is the cyclotron frequency, and 
*& n l(x, [(f)]) are the eigen-wave functions of the electrons under the external magnetic field B 
for a definite boson field 4>(x). It is well-known that at <f>(x) = the eigen-equation (|T3|) has 
exact solutions, 

^ nl (x,[0])=^ nl (x) = a nl (z-^-) (z*-i-) e~ zz * (15) 



dz* J \ dz / 

where / = 0, 1, 2, .., L max = $/$o, 3> = BS is the total flux, $ = 2irhc/e is the flux quantum, 
z = (x — iy)/(2l B ), z* = (x + iy)/(2l B ), l B = (he/ (eB)) 1 / 2 is the magnetic length, and a n i a 
normalization constant. Therefore, the eigen-functionals ty n i(x, [0]) have the following exact 
form, 

V nl (x, M) = A nl i> nl (x)e Q ^ x >® (16) 

where A n \ is a normalization constant, and the phase field Q n i( x :0 satisfy (static) Eikonal- 
type equation, 

+ (^A + a ni ) ■ pj Q nl (x, + [ ^' 0]2 + &{x) = S n ,[0] (17) 

where a ni = (l/m)phiif) nl (x). 

Following the above same procedures, we have the exact expression of the ground state 
energy as the electron density field p(x), 

E g [p\=E + l - j d 2 xd 2 y ^(x-y)p(x)p(y)-2^(x)^(y)^e(- J E^[0])^^J (18) 

where E = J2 n i E n i[(j)})e n , and the boson field <j){x) is the functional of the electron 
density p(x), and determined by the equation, 

Y,e(-E nl [<t>\) (G nl (x) + J d'y^y) 5 -^^) = (19) 
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where G n i(x) = J d£\ip n i(x)\ 2 e 2Q ni( x >® / J d 2 x\ijj n i(x)\ 2 e 2Q ^ x '^ , and the self-energy can writ- 
ten as a simple form, £ n j[0] = / d 2 x(f)(x)G n i(x). The chemical potential is determined by 
the constraint condition, 

J20(-E nl [<P]) = N (20) 

nl 

where N is the total electron number. The equations (|TTD, (|T5D , flT5D and (^) can be used 
to exactly determine the ground state energy and the ground state electron density of the 
system by taking SE g [p]/Sp(x) = 0. However, the expression of the ground state energy 
(PI) has a great advantages comparing with the Laughlin's trial wave-functions: 1). it is a 
microscopic theory expression. 2). it shows that the odd- and even-denominator's fractional 
quantum Hall states have the same expression of the ground state energy. 3). it can exactly 
determine the ground state electron density. 4). it is very simple to study the fractional 
quantum Hall effects taken place in higher Landau levels (n > 1). 

C. Electron gas with transverse gauge interaction 

We consider a two-dimensional electron gas represented by the Hamiltonian with trans- 
verse gauge fields A(x) (V • A(x) = 0), 



H = ^(x 



— (p-qA) 2 - a 
2rn y ; P 



ip(x)+~ I d 2 yv[x - !j)f)U)f)(!j) (21 j 



It is well-known that this system has the following low energy behavior, (1) as v(x — y) — 0, 
it shows the non- Fermi liquid behavior; (2) as A(a;) = and v (x) ~ In |x|, it also shows the 
non-Fermi liquid behavior. Thus the system belongs to the non-Fermi liquid. Introducing 
the Lagrangian multiplier 4>(x), we can re- write it as, 

1 



H[<f>,p]=^(x) 



(p-gA) 2 -p + <j){x) 
2m 



-(f)(x)p(x)+^ J d 2 yv(x - y)p(x)p{y) (22) 

Then by solving the eigen-functional equation of the "propagator" operator M(x) = 
(p-#A) 2 /(2m) -p + 4>(x), 
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2 i-(p-sA) 2 - p + * k (x, [0]) = £ fe [0]* fc (x, [0]) (23) 

we can obtain the exact functional expression of the ground state energy. The solutions of 
this equation are that, 

E k [(j>,A] =e fc + S fe [0,A] (24) 
S fc [0,A]= jf dij d 2 x%(x,[m[<t>(x) - ±M2p-gA)]* k (x,[t<l>]) 

where e k = {hk) 2 /{2m) — p, and A k is a normalization constant. The phase fields Q k (x,£) 
satisfy the (static) Eikonal-type equation, 



\2m mm J 



2m 



+ i{9 f* )? - ■ A(x) + &(x) = S fc [0, A] (25) 

and are the functional of the boson field <f>(x) and the gauge fields A(x). In general, this 
equation cannot be analytically solved but a series solution. Using the eigen-functionals 
ty k (x, to represent the electron field operators ip(x) and i]A{x), we can obtain the exact 
expression of ground state energy as a functional of the boson field 4>(x), the gauge field 
A(x) and the electron density p(x), 

E[(P,A,p]= J d 2 x<H[(P,p]> 

= XI 6(- E k[<f>, A])E k [(j>, A] - j d 2 x4>{x)p{x) + l - j d 2 xd 2 yv(x-y)p(x)p(y) (26) 
k J J 

where the boson field <p(x) is determined by the condition 5E[(j>, A, p]/S(j)(x) = 0, which leads 
to the equation, 

where G fc (x) = f^e^ x ^ / J <Pxe^ x ^ , and Q fe (x,0 = Qffof) +iQ^(x,0- The 
gauge fields A(x) are determined by the conditions Ji{x) = 5E[<f>, A, p]/5A i {x) — 0, i — 1, 2, 
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which means that the current induced by the gauge fields is zero. The ground state energy 
E[<f), A, p] ( |2SD has the similar functional expression to that of two-dimensional electron gas 
without gauge fields (|8|). However, they have different phase field dependence, this can be 
clearly seen from their self-energy expressionsE^ [(f)] and Efc[</>, AJ, respectively. The gauge 

fields produce the term %{x, [^4>})A(x)-p9 k (x, [£</>]) = \^ k (x, [£0])| 2 A(x)-(/ik+pQ fc (a;, £)) 
in the self-energy Efc[</>, A], thus in this case the imaginary part of the phase fields may 
contribute to the ground state energy. We would like to point that due to the phase fields 
Qk(%, are complex, to keep the eigen- values Ek[(j), A] being real, we must have an additional 
constraint condition on the gauge fields A(x)-pQ k = 0. Thus we have three equations 
to determine two independent transverse gauge fields. If these equations are independent, 
we only have the solution A(x) = 0, i.e., the pure transverse gauge fields do not contribute 
to the ground state energy. If they are not independent each other, we may have non-zero 
solution, and the transverse gauge fields contribute to the ground state energy. 

D. spin coupling systems 

Spin coupling systems are strongly correlated systems and generally represented by the 
lattice Hamiltonians. For simplicity, we only consider a D-dimensional spin-1/2 Heisenberg 
model, 

H = J^Si-Sj (28) 

<ij> 

where < ij > indicates the summation over the nearest neighbor sites. This spin-spin 
interaction represents the four-fermion/boson interaction, and is difficult to exactly treat it 
by usual perturbation theory. Here we only give a way for treating this kind of systems, 
which may give more useful and accurate informations than previous perturbation methods. 

In general, the spin operator (spin-1/2) Si can be represented by the fermion operators 
Si = \TioP^Jifi, and the Hamiltonian can be re-written as |T7| -§ ^2 <ij>)Cr)(r > fLfj*f] a 'fi*' 
or J2<ij>(fi^fji ~ fiifji)(fjlfa ~ fjlfil)) tllus we can write tne Hamiltonian as another 
form, 
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# = -^E^ + J E 5 iA> 

<ij> <ij> 



+J £ g . ^ + A * ■ & - \fL?af>fifi) (29) 



<ij> 



where xh — fiafjai an d = fjifi\~ fftfoi- The Lagrangian multiplier Aj take the equation 



Si = ^fJaCrapfifi as the constraint condition adding to the Hamiltonian (p9|), and the fermion 



operators f, ia and the spin operator Si are independent each other. In order to treat the 
four- fermion interactions in the Hamiltonian (p5|), we can further introduce the Lagrangian 
multipliers and Ay to decouple these interactions, and have the following form, 

H = J 2 E + J E 5 1A 

<ij> <ij> 

% + E A * • (S - \jl^^h) (30) 

<r/> i 

+ E M*c& - + h - c -] + E A ^ - S^t + Mn) + 

<ij> <ij> 

which is completely equivalent to the Hamiltonian (^8|). In fact, if we integrate out the 



boson fields Xij an d Dy, the Hamiltonian (|30f) becomes usual one []17|1 , and the Lagrangian 
multipliers and Ay become usual Hubbard- Stratonovich fields. 

Now the Hamiltonian (|30"D only includes the quadratic terms of the fermion operators, 
and in principle, we can exactly calculate the ground state energy. However, in the lattice 
sites, it is hard to obtain the eigen-functionals of the fermions and the ground state energy 
by present method, we need to mapping this lattice Hamiltonian into a continuous one. To 
this end, we take the following approximations, 

T* . I rp . ry . I ry . 

< Xij > = axi^ 1 ), < D i:j >= aD(^p) 



ry I ry ry I ry 

I)'.. — fit 1 ^ 3 \ \ \r > 



Ay = A(±i^i) (31) 



where a is the lattice constant, and <f)'( x ' ^ 3 ) = 7 + 4>( X *~^ X3 ), and obtain the continuous 
Hamiltonian, 
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H = J2e(k)fl k f ak - J d D x[<j>(x) + <t>*{x)]fl{x)Ux) + -J d D x\(x) ■ fl(x)a a0 f (x) 

+ J d D x{2A(x)ft(x)f l (x) - 2A*(x)fj(x)fl(x) + X(x)(fl(x)f a (x) - 1)} (32) 
+ J d D x{a.J\D(x)\ 2 - ^-\x{x)\ 2 } + J d D x[(p'(x)x(x) + A(x)D(x) + h.c] + H s 

where e(k) = -2^yJ2i C0S ( a k), H s = JJ2 <ij> S i -S j + J2 i Xi-Si, f ia = y/af a (xi), and X(x) is 
a Lagrangian multiplier which adds the constraint condition at each lattice site J} a fia = 1. 
The parameter 7 can be approximately determined by usual mean field theory as taking 
4>(x) = 0, where the parameter A(x) = Aq can also be determined by its self-consistent 
equation. In the continuous limit, the Hamiltonian Hg is reduced into one of usual non- 
linear a model under an "external magnetic field" X(x). 

We can further simplify the Hamiltonian (^) by taking e(k) = — /i, where m ~ I/7 
and \x is the chemical potential of the system. In general, due to A 7^ 0, there appears a 
gap in the excitation spectrum of the fermions, the dynamic term of the vector fields X(x) 
produced by the fermions is very small and can be neglected in the low energy (temperature) 
region ksT <C Ao, therefore, we can simply take A(x) = 0. We now briefly consider the 
influence of this approximation on the spin field. As J > 0, the system has a long-range 
or short-range antiferromagnetic correlation, thus the spin field operator can be written as 
Si = (— iySQ(xi) + aL(xi), where O(xj) • L(x«) = and |0(xj)| 2 = 1. In the continuous 
limit the coupling term becomes Yli^i ' &i ~ > ^ + ^/ a ) • O(g), thus only the 

fermion-hole excitations with large momentum (~ rr/a) contribute to this coupling term. 
However, due to Ao 7^ 0, these excitations are very small and can be safely neglected in the 
low energy limit. Therefore, the Hamiltonian (|32"D can be approximately separated into two 
parts, H = Hf + Hs, where Hf mainly represents the high energy (~ Ao) behavior, and Hs 
mainly represents the low energy (<^C Ao) behavior of the system. Here we only consider the 
contribution of the Hamiltonian Hf to the ground state energy, where the Hamiltonian Hf 
reads, 

H = J d D xfl(x){^ + X(x) - [0(x) + <F(x)]}f 9 (z) 
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+2 J d D x[A(x)f 1 (x)f l (x) - A*(x)fj(x)fl(x)} 

+ / d D x[—mx)\*-±\A(x)\*-\(x)] 
J aJ aJ 



(33) 



Due to appearing the pairing parameter A(x) = —aJ < f^(x)f^(x) >, the "propagator" 
operator M(x) of the fermions is a 2 x 2 matrix, and can be written as M(x) = M + $(x), 
where 

J - // -^-Afi 

M 



-2An 



2m 



( 
\ 



\ 



(34) 



\{x) - [(j)(x) + (f)*(x)] -2A (x) 
-2A(x) -X(x) + [<j)(x) + (f)*(x)} J 

where fx — /i — A , A(x) = A(x) — A , A(x) = \(x) — A , and the pairing parameter A 
can be determined by its self-consistent equation as taking <ft(x) = A(x) = \(x) = 0. The 
eigen-functional equation of the "propagator" operator M(x) reads, 



m(x)^\x, m) = E^[^:>(x, [$]) 



('<) r 



(35) 



where the eigen- values E^[$\ = ±E k + S^f^], E k = >Je 2 (k) + 4A§, and the self-energy 
= £ d£ / cPx^f > (x, [Z$\)$(x)¥f\x, [£$]). The eigen-functionals ^(x, [£$]) 



^(±)i 



can be generally written as, 



A (+) 



A 



(-) 



v k e c 



—vue Q< k +)( - x & 



„jkx 



/ 



■(-), 



„ikx 



(36) 



where w fc = (l/ v / 2)(l+£(A;)/ J B fc ) 1/2 , v k = (l/V2)(l-e(k)/E k )^ 2 , and are the normaliza- 
tion constants. To calculate the ground state energy, we only consider the eigen-functionals 
ty k \x, [£$]) which represent the states occupied by the fermions at zero temperature, and 



have the differential equations of the phase fields Q k \x,£) and Q k J (x,£), 
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where the pairing parameter A(x) is determined by its self-consistent equation A(x) = 
-aJ < f\(x)f\(x) >= ^ fc 0(-4 _) [$])^^e^" ) (^ +Q ^ ) *^. The functional expres- 
sion of the ground state energy can be written as a simple form, 

E[X,<j ) ,A] = J20(-El ) m)El ) m + I d D x[^\<f>'(x)\ 2 --L\A(x)\ 2 -\(x)] (38) 

k J 

If we can exactly solve the differential equations of the phase fields (|3~7D , we can exactly 
calculate the self-energy Sjj, '[<&] by its definition, then we obtain the ground state energy 
contributed by the fermions by taking the constraint conditions, 5 E[X,4>,A]/SX(x) = 5 
E[X, 0, A]/5(/)(x) = 5 E[X, 0, A]/5A(x) = 0. 

E. One-band Hubbard model 

We now consider the one-band Hubbard model on a square lattice, and give the exact 
expression of the ground state energy, which is very similar to that of the two-dimensional 
electron gas in the continuous coordinates. The Hamiltonian of the Hubbard model on the 
square lattice reads, 

H = -t ^a + A A) - v^a + u Y ( 39 ) 

<ij>,a icr i 

where < ij > indicates the summation over the nearest neighbor sites, and U is the on-site 
Coulomb interaction. To decouple the four-fermion interaction, we introduce the Lagrangian 
multipliers (fii a which take n ia = cf^c^ as the constraints to the system, and re-write the 
Hamiltonian ([39]) as, 

H[(f), n} = -t + c) a c ia ) - Y $A ~ Yl M n ™ ~ ^A) + u Y n ^ n ^ ( 4 °) 

<ij>,a icr icr i 

On the lattice site coordinates, the "propagator" operator of the electrons now can be written 
as a N x N matrix, 
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M ija ) = [M^ +<f> ia 5 i:j 



(41) 



where M„ 



(0) 



—t(jij + jji) — fi5ij, jij = 1 for j = i ± 1, and 7y = for other j. The 



eigen-functional equation of M a reads, 



(42) 



where the eigen-values E kfT [4>] = e k + Efc CT [0], e k = —2t(cos(ak x ) + cos(ak y )) — /z, the self- 
energy E fcCT [0] = Jq ^^[f^]^^^], and CT = (0i<j<%) is a N x N matrix, in which only 
the diagonal terms are non-zero. In general, the eigen-functionals Nt^^] can be written as 
a simple form, 



/ 



=1,5) \ 



p ik-Xjv p Qka{XN,£.) I 



(43) 



where A k is the normalization constant. With the equation ([42]), the phase fields Q ka (xi,^) 
satisfy the following equation, 



Ik 



I _ e -QUxi,0+QUx i+ i,t)} + 7 *(! _ e -Qk*(*i,0+Qto(*i-i,0) + ^ = s fc(J [0] (44) 



where j k = t(e mkx + e mhy ), and the self-energy £foj[0] is self-consistently determined by its 
definition. The electron operators Ci(j Cclll be represented by the eigen-functionals \l/fco-[</>] 
which are the eigen-functional wave functions of the electrons for a definite boson field 



tJ2 a 



(45) 



then we obtain the exact expression of the ground state energy as a functional of the boson 
fields 4>i a and the electron density fields , 



E[<j>,n] = < H[(p,n) > 



(46) 



where the chemical potential of the system is determined by the constraint condition, 



Y J 0(-EkM) = N e (47) 

where N e is the total electron number. By taking SE[(j), n\/54> io - = 8E[(f),n\/5n icT = 0, we 
can exactly calculate the ground state energy and the ground state electron density of the 
one-band Hubbard model by using the equations ([44]), ([46]) and ([47]). We hope that our 
present results can provide a powerful and exact method to study the strongly correlated 
electron systems represented by the Hubbard model, and can give some important criterions 
by numerical calculations for analytical calculations. 

F. Boson systems 

For a boson system, such as the liquid 4 He, at zero temperature it has Bose-Einstein 
condensation, and the bosons only occupy the state of the momentum k = 0, thus the 
equation (|i~0]) is trivial, and the ground state energy can be written as, 

E g [p\ = \J d D xd D y[v(x - y)p(x)p(y) - 25(x - y)(f>{x)p{x)} (48) 

where we have taken /i = S o [0]. The equation (||) reduces into, 

G (x) + j d D yHy) 6 -^y = ^P{x) (49) 

It is noted that the last term in (^) is the contributions of kinetic energy of the bosons, 
and in general it may be non-zero for interacting boson systems. We can also give the exact 
functional expression of the ground state energy of some trapped boson systems, where 
the procedures are very similar that we study the two-dimensional electron gas under the 
external magnetic field: (1). we first solve the eigen- functional equation of the "propagator" 
operator M(x) of the bosons under some external field and/or constraints, and obtain its 
eigen- values and the differential equations of the phase fields; (2). using the eigen-functionals 
to represent the boson field operators, we obtain the exact expression of the ground state 
energy as a functional of the Lagrangian multiplier boson field <fi{x) and the boson density 
field p(x). 
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III. Exact functional expression of the action 



In this section, we mainly focus on to study two kinds of systems, one is a D-dimensional 
electron gas with transverse gauge fields, and another one is the one-band Hubbard model on 
a square lattice. In these two cases, we shall give the exact functional expression of the action 
and the eigen-functional wave functions of the electrons. Under linearization approximation 
(only keeping the linear terms in solving the differential equation of the phase fields), we 
give the effective action which is equivalent to that obtained by usual RPA method. With 
them, we can calculate the excitation spectrum and the variety of correlation functions of 
the systems, and study their low energy behavior. 



A. Electron gas with transverse gauge field 



In general, we consider the system with the Hamiltonian (omitting spin label a of the 
electron operators C&) 

where A(x) (V • A(x) = 0) are the transverse gauge fields. The last term represents 
usual electron's Coulomb interaction (four-fermion interaction). Introducing the Hubbard- 
Stratonovich (HS) field <p(x, t) to decouple this four-fermion interaction term, and using the 
standard path integral method |[L8||, we have the action of the system, 



S[<f>,A] = dt d d x\^{x,t) 



at 2m 



+ 2^S^H.-fiM^) (51) 

where ^>(x,t) is the electron field, and only quadratic term of ty(x,t) appears in the action. 
Therefore, the integration of the electron field becomes a standard Gaussian form. It is 
worthily noted that after introducing the boson field <f)(x, t), the Hilbert space of the system 
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is enlarged, and the unphysical electron density field p(x,t) ^ ^(x,t)^f(x,t) appears in 
the action [H]. To erase the unphysical electron density field, we must take the functional 
average over the boson field (f>(x, t) in calculating a variety of correlation functions and 
response functions of the system to external fields. 

After integrating out t), the contribution of "kinetic" energy of the electrons to the 
action is —iTr ln(M), where M(x, t) = ihdt +fi — (p — gA) 2 / (2m) + 0(x, t) is the propagator 
operator of the electrons. Using the formula, 

Tr In (m) = Tr In (ihd t + fi-0 2 /(2m)) + [ d£ f dtd D x(p(x, t)G(x, t; x , t, [£</>]){ t ^ t 

\ ' Jo J x'—>x 

(52) 

where M(x,t)G(x,t; x ,t' , [(f)]) = <5(x — x.')5(t — t'), and neglecting the constant term, we 
have the following expression of the action, 

S[(j>,A] = -i [ d£ [ dtd D x(f>(x,t)G(x,t;x ,t , | t ,^ t 
+ 



o 



2TIP ^ v(q) 



The eigen-functional equation of the propagator operator M(x, t) reads, 

^ku(x, t, [(f)]) = E kuj [(f)]^ kw (x, t, [(f)]) (54) 



+ ^ ~ 2m ^ ~ gA ^ + ^ 



where the eigen-values Ek w [4>] = foui— £(&)+£& [0], e(k) = (Kk) 2 / (2m)— /x, and the self-energy 
S fc [0] = J d£ J dtd D x^l(x, [€4>])[<l>(x) + ^A-(2p-gA)]^f k (x, [£(/)]) is a regular function, and 
independent of u (see below equation |5(J). Using the orthogonality and completeness of the 
eigen-functionals *& ku> (x,t, [0]), the Green's functionals G(x,t;x ,t , [(f)]) can be represented 
by these eigen-functionals, 



G(x,t;x',t, [£</>]) = iriM*>»>(x,t, M)^Ux,t, [&]) (55) 

The eigen-functionals ^/ kw (x, t, [(f)]) are the eigen-functional wave functions of the electrons 
for a definite HS field (f>(x,t), and have the following exact expressions, 
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where is the normalization constant, and T — > oo is the time length of the sys- 
tem (for finite T, the ^kw(x, t, [£4>]) satisfy the boundary condition ^^(x, t + T, [£</>]) = 
— ^kio( x i t, These eigen-functionals are composed of two parts, one represents the 

free electrons, and another one represents the correlation of the electrons produced by the 
electron interaction. Thus we can formally write 

*Ux,t, [£</>]) = VUM)^e^«e- ffi *[* (57) 

where ^^{x^t) are the wave-functions of the free electrons. The phase fields Qk{x,t,^) 
satisfy the usual Eikonal-type equation with the condition Qk(x, t, £) = as <p(x, t) — and 
A(x,t) = 0, 

*^ar - o ( — k A ) ' P W x ' 



<9t 2m mm J 2m 

-^ + ^k.A + ^(x,t) = (58) 
2m m 

which can be exactly solved by a series expansion of the HS field 4>(x, t) and the transverse 
gauge fields A(x,t) and/or by computer calculations fl9|]2T)|| . It is noted that the <f>(x,t) 
and A(x,t) dependence of the eigen-functionals *$?(x, t, [£</>]) is completely determined by 
the phase fields Qk(x, t,£). 



Substituting equation (|56|) into equation (|55D , the Green's functionals can be written as 



G(x,t;x \t ',[&]) = J^0(£ fc [0] - £ (A;))|A fe |V k -( x - x ')-^( i - i ')e^^ + ^^' i '^ (59) 

k 

In order to calculate the Green's functionals G(x,t; x ,t , [£<f>] )| t /_ t , we take the following 

x'^x 

regular procedure, 

G(x,t;x',t',[Z4>])\ t ,^ = - lim (G(x,t;x-r],t,^(f)}) + G(x,t;x + r],t,[^})) (60) 
then we obtain the following exact functional expression of the action, 
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S[(f>, A ) = ^J dtdDx i F ( x > x-r},t)+ F(x, x + r), t)]„^ 0+ 

+ 22^E^(-9" n )^ n ) (61) 

f x ^\ e ih-(K-x.') e Q k (xM)+QU^'M) 



F(x,x>,t) = T[ dtJ2 9 &M-<k))^f 



dtd® xeQ k ( x ' t '^ + Q*k( x,t & 

It becomes very clear that this kind of problems end in to solve the differential equation fl58|) 
of the phase fields Qk(x,t,^), which can be solved by a series expansion of the boson field 
4>(x, t) and the transverse gauge fields A(x, t). With this action and the eigen-functionals, we 
can calculate a variety of correlation functions, such as the single electron Green's function, 
by taking the functional average on the boson field (f>(x, t) and the transverse gauge fields 
A(x,t), 



G(x — x' ,t — t') — % < G(x, t; x , t , [0]) ><£,a 

• / IW II DAG(x, t; x',t', [<j)])e iS ^ 



(62) 



/ 11 D<j> EI DAe lS ^ 

In general, it is difficult to exactly calculate the action, in order to study the low energy be- 
havior of the system, we can approximately calculate it. After taking some approximations, 
we can obtain the following effective action (omitting constant terms), 

qU 41 1 <l>(-q,-n)<Kq,n) 

q ,n W 

+~ J dtd D x<j)(x, t) [Ftix, t, S) + F 2 (x, £)]^ (63) 

where Fxfatrf) = -(1/(2tt) d ) / d D k6(-s(k)) sin(k • 77)7/ • VQ{{x, t), F 2 (x,t) = 
(2/(2tt d )) / d D k6(-e(k))Q%(x,t). We have taken 0(£ fc [0] - e(k)) = 9(—e(k)) and 
e Q k (x,w+Qt(x',t,o + Q k ( x ,t,£) + Q* k {x',t,£), and written the phase fields as Q k (x,t,£ = 
1) = Qk(x,t) = Qk (x,t) + iQl(x,t). It is worthily noted that the contribution of the "ki- 
netic" energy of the electrons to the action is composed of two parts, one is from the real 
phase field Q k {x, t), and another one is from the imaginary phase field Q{(x, t). Due to the 
Fermi surface structure of the system, the momentum integral in Fi(x, t, rj) can be written as 
% jd D k = So-! j d\k\\k\ D ~ l £ d6{sm6) D ~ 2 , where S D = 2ix D l 2 /T{D /2). As D > 2, the 
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integration of sin(k • 77) is regular, thus as 77 — > the function Fi(x, t, 77) = 0, the imaginary 
phase field Q{(x, t) does not contribute to the action. Only at D = 1, it has contribution to 
the action, where the real phase field Q§(x,t) is zero (see below). This property is indepen- 
dent of the electron interaction, it is completely determined by the Fermi surface structure 
of the system. 

For a ID interacting electron gas (taking A(x,t) = 0), the electron energy spectrum 
near its Fermi level ±kp can be written as, e(k) = ±v^k, where Vp is the Fermi ve- 
locity. The branch e(k) = vpk represents the right-moving electrons, and the branch 
e(k) = —vpk the left-moving electrons. In general, the electron interaction term reads, 
(1/-^) J2 q v ( ( l)PR( ( l)PL(— ( l), where Pr(l)(q) are the right- and left-moving electron densities, 
respectively, and v(q) = V, a constant. To decouple this four-fermion interaction, we in- 
troduce two HS fields <f>R(L)(x,t), an d have two sets of wave-functions of the right- and 
left-moving electrons for the definite <Pr(l)(x, t), respectively, 

^r( L )Ux, t, [&]) = I — J e Qn (L) (*M) e ik*-i{*-?; R(L) m (64) 

where T,r(l)[4>] is a regular quantity, and independent of k and uj. The phase fields 
QR(L)(%)t)£) are independent of k, and satisfy the simplified Eikonal equation, 
f d d\ 

\di ±iVF dx) ^*wOM'£) ~ €<f>R(L){x,t) = (65) 
These linear differential equations can be easily solved, and the phase fields QR(L)(x,t,^ = 
1) = Qn(L)(x,t) are imaginary because the HS fields <pR(L)( x ,t) are real. The action ( BTf ) 
becomes the following simple form . 



TL 

q,n 



9 \a 12 , 1 9 



4irq-n^ m ' 4vrg + fi'' r " 1 ' V 
where 4>R{L) — 0fl(L)(<2S It is worthily noted that the imaginary phase fields Qjjm (ar, t) 
not only determine the electron correlation, but also contribute to the action. This is 
qualitatively different from that in 2D and 3D electron gases, where the imaginary part 
of the phase field Qk(x,t) does not contribute to the action due to their Fermi surface 
structures. 
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(66) 



As D > 2, for simplicity, we can solve the Eikonal equation by neglecting the quadratic 
term (VQk) 2 - This approximation is reasonable because for the long-range Coulomb in- 
teraction, only the states near the Fermi surface with momentum q < q c (q c <C kp) are 
important in the low energy regime, and for the smooth function Qk(x,t), this quadratic 
term is proportional to {qc/kp) 2 ~ 0. Under this approximation, we can obtain the effective 
action, 



Seff. ~ 2TLD J2 



) |0(g,fi)| 2 + n i:; -(g,fi)A 4 (-g,-fi)^(g,fi) 



(67) 



where x(<7> ^) is usual Lindhard function [F23J, and Hij(q, fi) = (ryfi / 'q+XQ 2 )(^ij ~ QiQj I 'l 2 ) the 
propagator of the transverse gauge fields Ai(q, fl) produced by the electron- hole excitations 
|l7fl . In fact, the above approximation is equivalent to usual random-phase approximation 
(RPA). However, at present framework, it gives more useful and important informations 
than usual RPA method, because we have the well-defined phase fields, and can use their 
imaginary part to determine the electron correlation. 



B. One-band Hubbard model on a square lattice 



In the above subsection, we have considered the D-dimensional electron gas with the 
transverse gauge fields in the continuous coordinate space, now for simplicity we consider 
an one-band Hubbard model on a square lattice, but this method can be easily extended to 
other lattice models. The one-band Hubbard model on a square lattice is presented by the 
Hamiltonian, 

H = ~t ^ (CiaCjcr + c] a C ia ) -fiYl ^io-Cia + U ^ (68) 
<ij>a icr i 

Introducing the Lagrangian multiplier fields 0j ff (t), we have the following action, 

S[<l>,n]= Yl J dtc l(t)Mija{t)c ja {t) -J2 J dt[n ia {t)(f> i(T {t) + itfn^n^t)] (69) 

<ij>a ia 
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where Mij a (t) = Sij[id t + /i + 4>ia{t)] + ~t{lij + Iji) is a N x N matrix, where 7^ = 1 for 
j — i ± 1, and 7^ = for other j. After integrating out the electron fields c ia {t), we obtain 
the following action, 

S[<t>,n] = -iTr\n (M ija (t)) -J^J dt[n ia (t)(f) ia (t) + hjn^n^t)] (70) 

icr 

Using the formula, 

-iTr\n{M ijc {t)) = -iTr In (m^) - i^ i(J ft ^ / dt^^G^, t; Xj , t', [&]) \ t ^ t , where 
M-^ = 5ij(id t + /j,) + t(^ij + 7jj), and M i j cr (t)G cr (xj,t; x h t', [(f)]) = 5u5(t — t'), this action can 
be re- written as (neglecting the constant term), 

JO J xj—>Xi 

la 

~J2f dt[n ia (t)(f> ia (t) + If/n^n^)] (71) 
The eigen-functional equation of the matrix operator Mij a (t) reads, 



M a (t)V kuia (t, [(f)]) = E kuia [4>}y kuia (t, [(f)]) 



(72) 



where the eigen- values £^[^1] = a; — + £fc CT [0], £fc = — 2t[cos(ak x ) + cos(ak y )] — yU, and 
Zka[<t>] = Jo d£ J dWl wa (t, [4>])4>a{t)^ kuxT {t, [0]), where <j> a (t) = (<M*)5tf) is a iV x N matrix. 
The eigen-functionals ^ k(Jja {t, [(f)]) have the following general expression, 



1 \ 1/2 



^ e* k ' Xl e < ^ fcCT< - :El '*^^ ^ 



V 



,-i(u>+E fc<T [<0)i 



(73) 



/ 



e jk-xjv e Q fe(T (a;]v,i,5) 

where A fc is a normalization constant, and the phase fields Qka( x i,t,0 satisfy the following 
differential equation, 



_^y*(l _ e -Qfc<r( a; l ^.5)+Qfc<T( a; i-1^.5)^ — Q 



(74) 



where 7^ = t(e mkx +e taky ). Using the orthogonality and completeness of the eigen-functionals 
^kujait, [£</>]), the Green's functionals can be written as, 
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G ff ( Xi , t; Xj , t', M) = ^J2 6 ^ k ^ - 



■k 



2 ik-(x i -x J )-ie(A;)(t-t') 



k 



(75) 



then the action of the system reads, 




dt{n ia (t)<f) i(T (t) + -Unn(t)nn(t)] 



(76) 



icr 



F a {x i ,x j ,t)=T d£Y, e (Pk[<f>]-e(k)) 




It is clear that this action is very similar to that of the electron gas in the continuous coordi- 
nate space, the difference between them is that the phase fields Qkcr{%i,t,£) and Qk{x,t,l;) 
satisfy different differential equations ([7|) and (0), respectively. With equations Q73D, ( |7ip 
and (|76"1), we can completely determine the low energy behavior of the one-band Hubbard 
model on the square lattice, at least we can do that by numerical calculations, because we 
believe that with the help of these equations we may need much less computer time than 
previous methods for the large lattice site numbers. 

Now there is a common consensus that the physical properties of the high Tc cuprate 
superconductivity is determined by the 2D strongly correlated electron gas in their copper- 
oxide plane (s) which can be approximately represented by the one-band Hubbard model 
[/J. If this is true, our present results can provide enough informations to determine the 
physical properties of the high Tc cuprate superconductivity, because using equation ([16]) . 
we can exactly determine their ground state energy and ground state electron density field, 
and using equations ( |75| ) and ([7!]), we can calculate their variety of correlation functions 
and response functions to external fields. 



IV. Ground state wave function of the systems 



In principle, we can calculate the ground state wave function of any systems by using 
the eigen-functionals. For boson systems, at zero temperature they have the Bose-Einstein 
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condensation, and the bosons occupy the state with the momentum k = 0, and have the 
eigen-functional wave function ^?o(x, [</>])■ However, for fermion systems, the phase fields 
may depend on the momentum k of the fermions, it becomes very difficult to calculate their 
ground state wave functions. In general, the ground state wave function of the fermion 
system can be written as, 

*(ar 1 ,x 2 , ...,x N ) =< \\^ ki ( x j, [</>])\\ ></, (77) 

where k^ = kp is the Fermi momentum, and the functional average over the HS boson field 
4>(x) can be done by taking QkX x j) = Qhi x j^) an d 4>{ x i) = 4>(xi,t). Here we only consider 
the boson systems and the two-dimensional electron gas under an external magnetic field 
there appears the fractional quantum Hall effect for enough strong magnetic field. 



A. Boson systems 



For boson systems, due to the condensation of the bosons, we have the following ground 
state wave-function, 

*(x x ,x 2 , ...,x N ) = (7^72) <^(x 1 ,x 2 ,...,x N ,[(p])> ( f > 
^(xi,x 2 ,...,x N ,[(f>])= e ££i 00(^=1) (78) 

where < ... >^ means the functional average over the boson field <j>(x). The functional 



ty(xi, X2, xn, [(f)]) is very similar to the generalized London wave-function [24], where 
f(xi) = exp{Qj^(a^, £ = 1)} and S(xi) = Qb(xi,£, = 1). To calculate this functional average, 
we need knowing the action of the systems. Without external fields and constraints, the 
phase fields Qo(x,t,£) satisfy the simple Eikonal equation with the condition Q Q (x,t,£) = 
as <f>(x, t) = 0, 

%n dt-^) g ° (X ' *' 2^ + ^ t] = ° (?9) 

If we only keep the linear terms, we have the following simple solution, 
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E 



Ml — e q 



(80) 



where e g = h 2 q 2 /(2m). Substituting it into fl63|) , we obtain the following effective action 



with k = 



5*0 



(M)) 2 -e2 2u(g) 



iVe 9 1 




(81) 



It is worthily noted that with this action we cannot obtain the correct excitation spectrum 
of the systems, because it is obtained by using the phase fields Qo(x,t,£). In order to get 
the correct excitation spectrum of the bosons, we need knowing the action obtained by the 
phase fields Qk(x,t,£). However, we can approximately use the action So[4>] to calculate the 
functional average over the boson field <fi(x,t) in equation ( |78| ) by taking Qo(%i,£, = 1) = 
Qo(xi,t,£ = 1). The detail expression form of the ground state wave function depends on 
the interaction potential v(q), but it can be generally written as, 



where U(xi — Xj) is uniquely determined by the interaction potential v(xj — Xj). It is clearly 
seen that this ground state wave-function is uniquely determined by single effective potential 
function, and has the expression very similar to usual correlated basis functions [?] that are 
the type of wave-function most often employed in the study of the ground state properties 



The ground state wave function of this system cannot be written as a simple form even 
if for the lowest Landau level (n = 0), because the phase fields Q n i(x,£ = 1) depend 
on the quantum numbers n and I. In the lowest Landau level, if we approximately take 
Qoi{ x ii = 1) ~ Qoo(^), we can obtain the following expression of the ground state wave- 
function, 



x 2 , ...,x N ) ~ e 



Ei<j U{xi-Xj) 



(82) 



of 4 He. 



B. Two-dimensional electron gas under an external magnetic field 
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JV 

V( Xl ,x 2 ,...,x N ) = Ajlzi Hizi-z^tiM 2 <e Ef =1 QooK) >(j> (83) 

i=l i>j 

where A is a normalization constant, and the factor n^.-^i ~~ z j) guarantees the anti- 
commutation of the electrons. The last factor < exp{J^ =1 QooO^i)} ><t> is very similar to 
the ground state wave- function of the boson systems (0), and it is very clear that this factor 
is the contribution of the electron interactions. Using the procedure as above for the boson 
systems, we can obtain an approximate action by only keeping the linear terms in solving 
the differential equation of the phase field Qoo(^) ; then with this action we can calculate the 
functional average in (B3). Only for long range Coulomb interaction v(q) = e 2 /(4nq 2 ) can 
one have < 6xp{J^ i=1 Qoo(xi)} ></,~ Yli>j \ z i ~ z jV ■> where 7 is a dimensionless constant, 
but at present approximation, it is hard to determine the relation between 7 and the filling 
factor of the lowest Landau level. Thus this ground state wave function is different from 
the Laughlin's trial wave-functions of the odd-denominator fractional quantum Hall states. 
However, we believe that with the equations of the ground state energy and the ground state 
wave function ([18]) and (0), we can obtain more important informations of the fractional 
quantum Hall effects than that by the Laughlin's trial wave-functions, because they are 
directly derived from the microscopic theory. 

V. Unified description of strongly and weakly correlated electron gases 

In this section, we show that the strongly and weakly correlated electron gases can be 
unifiably represented under present theoretical framework. In fact, as shown in the section II 
and /or section III, the problems of the quantum many-particle systems end in to solve the 
differential equation of the phase fields Qk(x,t,£). Thus the phase fields Qk(x,t,£) are the 
key parameters hidden in the quantum many-particle systems, and completely determine 
their low energy behavior. It is natural that we can unifiably represent the quantum many- 
particle systems with the help of the phase fields Qk(x,t,£). 
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It is simple to prove that the ID interacting electron gas is a strongly correlated system 
even for very weak electron interaction V ~ 0. With the action (|66|), by simple calculation 
we can obtain the relations, 

< ML)Ux,t)^ R{L)k Jx',t, [0]) ><p ~ (lj e *(*-'> 

< e O«(z.)(x,<) e -0«(z.)(x',t) > ^ e -2a\n\x-x'\ ^ ^ _ ^ _^ ^ ^ 

where a ~ (1 — V/(2irhvF))/2 for V ~ 0, is the dimensionless coupling strength param- 
eter, the if>R(L)kw(x,t) are the wave-functions of the right (left)-moving free electrons, and 
< ■ ■ • >0 means taking functional average over the HS fields (pR(L)(x,t). The first equation 
presents the zero overlap between the eigen- functional wave functions ^ R(L)kuj(x,t, [(f)]) of 
the interaction electrons and the wave functions of the free electrons (V = 0) as L — > oo, 
thus the states of the interacting electron gas does not have one-to-one correspondence via 
adiabatic continuation with those of the free electron gas. The second equation presents 
the strong electron correlation, in which the electron correlation length is infinity even for 
weak electron interaction. Due to this strong electron correlation, the low energy excitation 
modes of the ID interacting electron gas are those collective excitation modes, such as the 
charge and spin density waves, and there are not well-defined quasi-particles (holes) near the 
two Fermi levels ±/cf. Thus the Tomonaga-Luttinger liquid theory is a strongly correlated 
theory, and is universal for ID interacting electron gases. 

For a 3D electron gas with long-range Coulomb interaction v(q) = e 2 /(Airq 2 ), by simple 
calculation, we have the relations, 

< ^(x,t)^l(x',t, [0]) ><t> ~ Z k e^-^ 

< e iQl(x,t) e -iQ{(x' ,t) > ^ ^ e z' k (x-x') ^ 

where Z k is finite, and z T k {x) is a smooth function. As |x| — > oo, the function z T k (x) goes 
to zero. The first equation presents that the eigen-functionals ^kuj(x,t, [0]) have large over- 
lap with the eigen-functions ^^(rr, t), and the second equation presents that there is only 
weak electron correlation even for long-range Coulomb interaction. Thus the fundamental 
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assumption of the Landau Fermi liquid theory is satisfied, and the Landau Fermi liquid 
theory is universal for 3D interacting electron gases. 

For a 2D electron gas with the long-range Coulomb interaction v(q) = e 2 /(4irq 2 ), the 
situation is different from that in the 3D electron gas. By simple calculation, we can obtain 
the relations, 

/9 



t fV + U1W , ~ I — I **-(*- x ') 



< ^(x,t)^( x ' 5 t ; [0]) > ^ i _ i ( 



1 

7 

< e iQU*,t) e -iQl(x',t) > ^ „ e -2/3M 9 c|x-x'|) ) | x _ x ,| ^ ^ (g6) 



where j3 = e 2 /(2(4ir) 2 ujp) is the dimensionless coupling strength parameter, and uj p is the 
plasma frequency. In fact, q c L — > oo, the first equation presents that the eigen-functionals 
^kuj(x,t, [4>}) has no (or infinitesimal) overlap with the eigen-f unctions i^ku{x, t), and the 
second equation presents the strong electron correlation. In this case, the fundamental 
assumption of the Landau Fermi liquid theory fails, and the system shows the non-Landau- 
Fermi liquid behavior. In the low energy limit, as Q > q —>■ 0, we have x(<? ; ^) ~ Q 2 /^ 2 , then 
we obtain the propagator of the HS field <fi(x,t) by the action (p7|), G^q, Q) = v(q)/(l — 
uj 2 /VL 2 ). Using this propagator to calculate the functional average over t), we obtain the 
results in (]5^|). These results are consistent with that of Ref. ^E] obtained by other method. 
The higher order terms neglected in solving equation (|58D cannot alter these results because 
in low energy limit they are proportional to (q c /kF) 2+n — ► 0, n = 0, 1, 2, and can be safely 
neglected. It is worthy noted that for 2D and 3D electron gases with long-range Coulomb 
interactions we have the same G^(q, Q) in which the screening effects are denoted by the 
factor 1/(1 — u 2 /Q 2 ). However, the reason they show different low energy behavior is that 
the integral J d d qe lcl ' x /q 2 appearing in the functional average of equations (^5|) and ([86]) is 
divergent (log-type) for d = 2, while it is convergent for d = 3. 

Even if the asymptotic behavior of the single electron Green's function is similar to 
that of the ID interacting electron gas (see, ( |8"1D and (|86|)), they originate from different 
physics. However, this difference can also be clearly seen from their effective actions of the 
density field. For ID interacting electron gases, this action can be easily obtained by usual 
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bosonization method [ I3|,|15 |, where the density field has well-defined propagator (density 
wave). For the 2D electron gas with long-range Coulomb interaction, the effective action 
(|67|) can be written as with the density field (taking A(x,t) = 0), 

There appears a gap in the excitation spectrum of the density field (plasmon excitations). 
It was shown [2£| that as e 2 3> 47r the system can be described as a Landau Fermi liquid 
formed by chargeless quasi-particles which has vanishing wavefunction overlap with the bare 
electrons in the system. For short range and/or weak electron interaction, the 2D interacting 
electron gases would show the Landau Fermi liquid behavior in the low energy region. 

Based on the above calculations, it becomes clear that if the non-Fermi liquid behavior of 
the electron gases is produced by the (long range) Coulomb interaction of the electrons, the 
correlation function of the imaginary phase fields Q{{x, t) must be a log-type. However, there 
exists another type correlation function of the imaginary phase fields Q{(x,t), which also 
induces the non-Fermi liquid behavior of the electron gases. Here we only consider a simple 
example of 2D electron gas with transverse gauge fields A(x,t). Solving the differential 
equation of the phase fields Qk(x,t) (|5S[) where only keeps the linear terms, and using the 
effective action (|67D, we can obtain single electron Green's function, 



G(x - x',t-t') =i < G(x, t; x', t', [</>]) > A 

= - \ 9(-e(k))e^< x - x ' ) - ie ^ t - i/) e p ^ x - x '^ (88) 
k 

where P£(x,t) comes from the contribution of the imaginary phase fields Q ! k (x,t), and the 
contribution from the real phase fields Qf(x,t) can be neglected. Near the Fermi surface 
k ~ fop, we have the relation, 

Pi F (0,t) ^ -t 1/3 ag 2 t 1/3 (89) 

where a = ^^ (x/l) 1 ^ 3 J dx(l — cos(x))x~ 4 / 3 . This result is basic same as that of Ref. 
27| , |28[] . Due to the power-law time dependence of the phase factor P^ F (0,t), the single 
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electron Green's function shows a singular low energy dependence which violates the quasi- 
particle excitation assumption of the Landau Fermi liquid. Thus, the system shows non- 
Fermi liquid behavior. In this case, the correlation function of the imaginary phase fields 
Q{(x,t) has the power-law form t a , a < 1, near the Fermi surface. For a 2D electron gas 
with long range Coulomb interaction and transverse gauge fields, the correlation function of 
the imaginary phase fields Q{(x,t) is composed of two parts, one is a log- type which comes 
from the contribution of the HS boson field <f>(x, t), and another one has the power-law form 
t a , a < 1, near the Fermi surface, which is from the contribution of the transverse gauge 
fields A(x,t). Thus this system shows the non- Fermi liquid behavior in the low energy 
region. 

VI. Discussion and conclusion 

Applying previous perturbation methods to study strongly correlated systems, we may 
meet many serious problems, because in these systems due to the strong correlation among 
fermions/bosons, there does not exist a small quantity to be used as a perturbation ex- 
pansion parameter. The most prominent character of the strongly correlated systems is the 
strong fermion/boson correlation, which completely controls their low energy behavior. This 
strong fermion/boson correlation is produced by the fermion/boson interactions, thus the 
key point of this kind problems is how to exactly or accurately treat the particle interac- 
tion (potential) terms. However, the traditional treatment of these systems is starting from 
a bare particle Green's function to perturbatively treat the particle interaction (potential) 
terms, i.e., exactly treating the kinetic energy term and perturbatively treating the potential 
terms. Explicitly, this treatment is unsuccessful for the strongly correlated systems. 

It is natural that one uses another way to treat the strongly correlated systems, i.e., 
exactly treating the interaction (potential) terms and exactly or perturbatively treating the 
kinetic energy term of the particles, such as usual bosonization method in treating ID in- 
teracting fermion systems. Only in this way can one effectively treat the strong particle 
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correlation. However, to realize this idea, we may meet two problems, one is how to de- 
couple four-particle interaction (potential) term, and another one is how to represent the 
particle correlation strength. The former one can be done by introducing Lagrangian mul- 
tiplier or Hubbard- Stratonovich boson fields, and the latter one can be naturally done by 
introducing phase fields in solving the eigen-functional equation of the propagator operator 
of the particles. It must be mentioned that after introducing the Lagrangian multiplier or 
Hubbard-Stratonovich boson fields, the Hilbert space of the system is enlarged, and the 
unphysical particle density field also appears in the Hamiltonian and/or action. To delete 
these unphysical states or leave the original Hilbert space intact, we can add a constraint 
condition to the functional ground state energy, such as 5E[(f), p]/8<j)(x) = 0, or take the 
functional average over the Lagrangian multiplier <f>(x,t). 

Based on the above idea, we have established an unified theory of the quantum many- 
particle systems, which is valid not only for weakly correlated fermion/boson systems, but 
also for strongly correlated fermion/boson systems. This unified theory of the quantum 
many-particle systems has three prominent characters, (1). it is completely founded on 
the basic principles of quantum mechanics. (2). it has well-defined phase fields used to 
represent the particle correlation strength. (3). it has an explicit and simple exact functional 
expression of the ground state energy, the action and the eigen-functional wave functions of 
the particles, and one can easily use them to effectively (or exactly) treat general quantum 
many-particle problems, and can obtain any accurate result as one hoped by taking high 
order corrections. In this unified theory, we have shown that the problems of the quantum 
many-particle systems end in solving the differential equation of the phase fields Qk(x, 
and/or Qk(x, t, £), and given the exact expression of the ground state energy and action, 
and the eigen-functional wave functions of the particles as the functional of the HS boson 
fields (and/or transverse gauge fields) and the particle density fields. With them, we can 
calculate a variety of correlation functions of the systems. It is surprising but natural that the 
physical properties of the many-particle systems is completely controlled by the differential 
equation of the phase fields Qk(x,t, £). The above results have been easily applied to the 
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systems under external fields and/or constraint conditions, such as trapped boson systems 
and 2D electron gas under the external magnetic field, and to the lattice models, such as 
the Heisenberg model and the one-band Hubbard model on the square lattice. 

The key points of this unified theory are that: (1). we introduce the Lagrangian multiplier 
(or HS) field <j)(x,t) which takes the particle density p(x,t) = ipi(x,t)ip(x,t) as a constraint 
condition, so that the four-particle interaction term is decoupled and the action can only 
have the quadratic form of the particle (fermion/boson) fields. (2). by introducing the 
phase field Qk(x, t,£), which is a functional of the boson field <j>(x,t) and determined by 
usual Eikonal-type equation, we can use it to completely represent the "kinetic" energy of 
the systems. Thus the problems of the quantum many-particle systems end in to solve the 
differential equation of the phase fields Q k (x,t,£). (3). the phase field Q k (x,t,£) is a key 
parameter hidden in the quantum many-particle systems, its imaginary part represents the 
particle correlation strength, and its real part only contributes to the ground state energy and 
action. (4). we are able to use the particle density field p(x) and p(x, t) to exactly represent 
the ground state energy and action of the quantum many-particle systems, respectively. 
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